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LOCALIZATION OF CUSTOMER PREMISES IN A LOCAL LOOP 
BASED ON REFLECTOMETRY MEASUREMENTS 

The present invention relates to a method to qualify a local loop 
from reflectometry measurements thereof, as described in the preamble of 
5 claim 1, and to an apparatus wherein the method is applied, as described in 
the preamble of claim 6. 

Such a method and apparatus wherein such a method is applied 
are already known in the art, e.g. from the European patent application 
number 01400832.0 entitled 'Method and Apparatus for Identification of an 
10 Access Network by means of 1-port Measurements', being filed on 2 nd April 
2001 in the name of Alcatel. Therein, a method to predict network 
parameters of a local loop from reflectometry measurements thereof is 
described. 

A local loop is a two-conductor transmission line that connects 
15 customer premises to a central office of a telecommunication network. 
Originally, the local loop has given the customer access to the Plain Old 
Telephony Service (POTS). This copper-wired access network has been de- 
facto the medium to deliver at the best cost new telecommunication services 
to residential users. Digital Subscriber Line (DSL) services exploits the 
20 frequency band above 4 kHz up to several MHz, which is not used by the 
POTS, to provide the customer with high speed internet access (up to a few 
Mbit/s). 

Nevertheless, the local loop was engineered for voice-band 
transmission without any prospect for data telecommunication services 

25 operating at higher frequency. The local loop may have several impairments 
for DSL services that do not affect the POTS, such as bridged taps (open 
ended line derivation spliced or coupled to the main line), mixed wire gauge, 
splices, etc. DSL services are mostly affected by reflections arising along the 
line, which stem from any non-uniformity in the impedance of the 

30 transmission line. Reflections interfere with the direct DSL signal and, at 



- 2 - 



those frequencies for which the interference is destructive, reduce the 
magnitude of the received signal. 

It is therefore necessary to evaluate a priori the line quality 
before implementing DSL services. An object of the known method is to 
5 predict the maximum bit rate achievable over a local loop - also referred to 
as theoretical channel capacity of the local loop - for a telecommunication 
service to be delivered over the local loop with the expected quality, namely 
to qualify a local loop for a telecommunication service. 

In a preliminary step of the known method, a physical model of 
10 the local loop is defined from the electromagnetic theory specifying how the 
voltage and current waves propagates throughout the loop. The local loop is 
modeled as comprising a plurality of line sections spliced or coupled to each 
other in a tree-like structure wherein: 

- the root of the tree is a line termination being wired to the central office, 
15 - a branch of the tree is a two-conductor uniform line section, 

- a node of the tree is a connection between a plurality of line sections, 

- a leave of the tree is a line termination, all being left open except one 
being wired to the customer premises. 

The root or the leaves of the tree are referred to as loop termination further 
20 in the text. 

A line section is modeled by the brace (z o ,e - ^) where Zo is the 

characteristic impedance of the line section, ^is the propagation constant of 
the line section, and I is the length of the line section. Zo and y can be 
expressed as a function of the frequency / and of a finite number of 
25 coefficients, the coefficients being defined from the cross-sectional geometric 
dimensions of the line and the material constants of the line. (z o ,e~^) 

specifies the way the voltage and current waves propagate along the line as 
well as their respective ratio. The boundary conditions, which the line 
voltage and the line current must comply with at a splice or at a line 
30 termination, yield the reflections and the transmission factor of the voltage 
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and current waves at that splice or at that line termination. Extending the 
reasoning to the entire loop yields a parametric model specifying how the 
voltage and current waves propagate throughout the loop, the above 
mentioned coefficients being the parameters thereof. The network 
5 parameters of the local loop, such as the scattering matrix or the ABCD 
matrix, are computed therefrom. 

In a further step of the known method, an initial electrical signal 
is generated and coupled to a loop termination, while the electrical 
waveforms at that coupling point are observed. The observed waveform is 

10 comprised of an amount of the initial signal and an amount of each signal, if 
any, that is reflected back from any non- uniformity in the impedance of the 
transmission line. This measurement procedure is known in the art as 
reflectometry, more specifically as time domain reflectometry or frequency 
domain reflectometry depending on whether the observed waveform is 

15 measured in the time domain or in the frequency domain respectively. The 
measured signal is called a reflectogram, a time domain reflectogram or a 
frequency domain reflectogram. 

At that stage of the qualification process, the local loop is 
considered as a one-port system, the loop termination at which the 

20 reflectometry measurements are performed being the reference port thereof. 
The loop termination being wired to the customer premises is loaded with a 
first impedance, given a definition thereof (e.g., an open-circuit) or, at least, 
an approximation in the considered frequency domain (e.g., by means of a 
RLC circuit). 

25 In a further step of the known method, the coefficient values of 

the parametric model are estimated. The estimation process includes the 
steps of: 

- assigning an initial value to some of the coefficients by making use of a 
quasi-realistic model of the local loop, 
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- computing from the parametric model the scattering parameter S u of the 
local loop, namely S™ del , the scattering parameter S u being defined as 
the ratio in the frequency domain of the backward travelling wave or 
reflected wave over the forward travelling wave or incident wave, 

5 - computing from the reflectometry measurements the scattering 
parameter S n of the local loop, namely s™ easuremettt , 

- quantifying a difference between the scattering parameters S™ del and 

^measurement by introducing R CQst function, 

- reducing the difference by adapting the coefficient values of the 
10 parametric model. 

The quasi-realistic model provides initial values to the coefficients 
of the parametric model for the computation of the scattering parameter 
S modei The initial values of the coefficients should preferably have the order 
of magnitude predicted by their formula so that the optimization of the cost 

15 function determines a absolute minimum and does not get trapped in a local 
minimum. The quasi-realistic model is a white box model with as much a- 
priori knowledge as possible, as opposed to a black box model. The quasi- 
realistic model of the local loop is based upon an assumed or an estimated 
loop topology, and possibly upon the type of transmission line used in the 

20 field as presumably known to the telecommunication operator. A method to 
estimate the topology of a local loop from reflectometry measurements 
thereof is already known in the art, e.g. from the European patent 
application number 02291902.1 entitled Interpretation System for 
Interpreting Reflectometry Information', being filed on 26 th July 2002 in the 

25 name of Alcatel. Therein, an expert system built upon a belief network and 
a rule-based system is described that determines the tree structure of the 
local loop and estimates the delay introduced by each line section. 

Once the coefficients have been estimated, it is possible to 
compute the transfer function between any loop termination, the local loop 

30 now being considered as a multi-port system. Given the so-computed 
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transfer function between the central office and the customer premises, 
given an estimated or measured Power Spectral Density (PSD) of the noise 
at the customer premises, it is possible to compute the theoretical channel 
capacity of the local loop. 
5 While the loop termination associated with the central office can 

be easily ascertained, guessing which one of the loop terminations is 
associated with the customer premises is not that obvious when there are 
more than one candidate. Without that piece of information, it is impossible 
to determine which one of the transfer functions must be used to compute 

10 the theoretical channel capacity of the local loop. 

An object of the present invention is to determine which one of the 
loop terminations the customer premises are associated with, given an 
assumed or an estimated topology of the local loop. 

According to the invention, this object is achieved by the method 

15 defined in claim 1 and by the apparatus defined in claim 6. A second 
reflectogram is carried out while the loop termination being wired to the 
customer premises is loaded with a second impedance distinct from the first 
impedance. Next, the first reflectogram is compared in the time domain 
with the second reflectogram in order to delineate among all the reflections 

20 those at the customer premises. Indeed, an impedance change at one loop 
termination affects only the reflections at that loop termination, all the 
other reflections remaining unchanged. The observed waveforms diverge 
one from the other whenever signals reflected back from the customer 
premises are present. The first divergence peak delineates the first order 

25 reflection. Therefrom, the round trip delay between the measurement side 
and the customer side is estimated. It is noteworthy that the magnitude of 
the divergence peaks decreases as the time increases, since the more delay 
the more attenuation. The skilled person will appreciate that other features 
in the outcome of the comparison step could be used to estimate the round 

30 trip delay with an acceptable accuracy. 
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Given the following pieces of information, as provided by the 
telecommunication operator and/or by the above mentioned expert system: 

- the tree structure of the local loop, 

- the delay introduced by each line section, 

5 it is possible to select, out of all the paths between the loop termination at 
which the measurements are performed and any other loop termination, the 
path which the accumulated delay matches the closest the measured round 
trip delay. 

Once the loop termination associated with the customer premises 

10 has been identified, the proper transfer function can be selected for the 
computation of the theoretical channel capacity between the central office 
and the customer premises. 

Preferably, such an impedance switching at the customer side 
should be remotely controllable or remotely observable at the measurement 

15 side, so as to keep the whole qualification process at one side. Various 
embodiments for switching remotely the impedance load at the customer 
premises can be thought off. The claim 2 defines such an embodiment if the 
POTS is already supplied to the customer over the local loop. The first 
reflectogram is executed while all the phones at the customer premises are 

20 on-hook (very large impedance acting as an open circuit), the phone line 
being idle, and the second reflectogram is executed while one phone at least 
is off-hook, the phone line being busy. A maintenance phone call to the 
customer may initiate the impedance switching. 

Another characterizing embodiment of the present method and of 

25 the present apparatus is defined in the claims 3 and 7. The reflectometry 
measurements are performed in the time domain to proceed directly with 
the comparison step without any further processing, provided the same 
signal has been used to excite the loop. 

Another characterizing embodiment of the present method and of 

30 the present apparatus is defined in the claims 4 and 8. The reflectometry 
measurements are performed in the frequency domain, next the 
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measurement results are converted back to the time domain, e.g. by means 
of the Inverse Fast Fourrier Transform (IFFT), before proceeding with the 
comparison step. 

Another characterizing embodiment of the present method and of 
5 the present apparatus is defined in the claims 5 and 9. An incident and a 
reflected wave are delineated from the observed waveform, e.g. by 
measuring both the line voltage and the reflected voltage wave. It is then 
possible to compute the time domain reflectogram of a reference signal, such 
as a direct pulse. In so doing, the excitation signals need not to be identical 
10 to proceed with the comparison step. 

Further characterizing embodiments of the present method and of 
the present apparatus are mentioned in the appended claims. 

It is to be noticed that the term 'comprising', used in the claims, 
should not be interpreted as being restricted to the means listed thereafter. 
15 Thus, the scope of the expression 'a device comprising means A and B' should 
not be limited to devices consisting only of components A and B. It means that 
with respect to the present invention, the only relevant components of the 
device are A and B. 

Similarly, it is to be noticed that the term 'coupled', also used in the 
20 claims, should not be interpreted as being restricted to direct connections only. 
Thus, the scope of the expression "a device A coupled to a device B" should not 
be limited to devices or systems wherein an output of device A is directly 
connected to an input of device B. It means that there exists a path between 
an output of A and an input of B which may be a path including other devices 
25 or means. 

The above and other objects and features of the invention will 
become more apparent and the invention itself will be best understood by 
referring to the following description of an embodiment taken in conjunction 
with the accompanying drawings wherein: 
30 - fig. 1 represents an apparatus according to the invention, 
- fig. 2 represents an instance of a local loop, 
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- fig. 3 represents a first reflectogram (inc. the initial excitation signal), 
executed while all the phones at the customer premises are on-hook, 

- fig. 4 represents a second reflectogram (inc. the initial excitation signal), 
executed while one phone at least at the customer premises is off-hook, 

5 - fig. 5 represents the comparison of the first reflectogram with the second 
reflectogram, achieved by subtracting the second reflectogram from the 
first reflectogram. 

A preferred embodiment of the present invention as depicted in 
fig. 1 comprises: 
10 - an input device IN, 

- an output device OUT, 

- a time domain reflectometry unit TDR and 

- a processing unit PROC. 

The input device IN and the output device OUT are coupled to the 
15 processing unit PROC. The time domain reflectometry unit TDR is coupled 
to the processing unit PROC via a data communication bus COM, such as a 
serial bus, and to the local loop LP. 

The apparatus APR is integrated as a handheld testing device 
that can be used by a technician in the field, e.g. at the central office next to 
20 a distribution frame. 

The input means comprises a keyboard and a computer mouse 
associated with a graphical user interface and are adapted to: 

- trigger the time domain reflectometry measurements, 

- enter or select from a drop -down menu any useful pieces of information 
25 such as characteristics of the local loop to be qualified (used in the 

parametric model), measurement settings, etc. 

The output device OUT comprises a display screen to provide a 
technician with the outcome of the qualification procedure, such as a yes/no 
answer to whether the loop capacity is greater or less than a pre-determined 
30 target figure. 



The time domain reflectometry unit TDR comprises the following 
functional units: 

- a generation unit GEN, to generate an electrical signal of a given shape 
and within a given spectrum range, 

- a measurement unit MEAS, to measure in the time domain an electrical 
signal, such as a voltage, including an analog to digital converter, 

- a coupling unit COUP, to couple the excitation signal from the generation 
unit GEN to a conductor pair and to provide the measurement unit 
MEAS with the appropriate measurement access points, 

- a filtering unit FILT, to filter the signal from the local loop LP, 
comprising: 

- a high-pass filter, to attenuate strongly the frequency band up to 4 
kHz, which may carry POTS signalling and voice traffic, 

- an low-pass anti-alias filter for the analog to digital conversion, 

- a low-noise amplifier LNA, to amplify the signal from the local loop LP, 

- a synchronization unit SYN, to trigger synchronously the generation unit 
GEN and the measurement unit MEAS. 

The processing unit comprises a Pentium III processor supplied 
by Intel Corp, as well as all the necessary hardware as known to the skilled 
person, such as random access memory, a hard disk, an interruption 
controller, a DMA controller, a graphic card, I/O peripherals, etc. 

The processing unit PROC runs software which carries out any of 
the method in accordance with the present invention, more specifically: 

- gather and average the measurement results from the time domain 
reflectometry unit TDR, 

- estimate from these measurement results the topology of the local loop, 
as already known to the art, 

- estimate from these measurement results the network parameters of the 
local loop, as already known to the art, 

- determine from these measurement results which one of the loop 
terminations the customer premises are associated with, 
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- qualify the local loop for a specific telecommunication service. 

The electrical signal used to excite the local loop is a block pulse 
signal. The signal type as well as the burst duration can be set via the input 
means IN. 

5 The measured electrical signals are the line voltage VIN in 

conjunction with the reflected voltage wave VIN-. The latter allows the 
discrimination of the incident and reflected waves within the line voltage in 
order to compute the scattering parameter s™ easurement . 

Upon request from the operator via the input device IN, the 

10 processing unit PROC triggers the time domain reflectometry unit TDR via 
the data communication bus COM for it to carry out a time domain 
reflectometry measurement. Thereupon, the synchronization unit SYN 
triggers both the generation unit GEN, to generate an excitation signal, and 
the measurement unit MEAS, to start measuring the resulting signals at 

15 the coupling point. A measurement is a set of digital samples, each sample 
being the encoded magnitude of the measured signal at a particular 
instance of time. The measurement results are delivered via the data 
communication bus COM back to the processing unit PROC for further 
processing. The measurement procedure is repeated a few times over a 

20 small period so as to compute the mean and variance value of each sample. 

In the following, the operation of the preferred embodiment will 
be described with reference to an instance of a local loop. However, the 
present invention is not limited thereto. 

The local loop LP as depicted in fig. 2 comprises a single 

25 transmission line TL with at some point along the line a bridged tap BT. 
The apparatus APR stands at the central office CO and is coupled to the 
termination TERM 1 of the local loop LP. The apparatus APR is modeled as 
a voltage source E serially connected to an impedance ZG and as a 
measurement device M measuring the line voltage VI and the reflected 

30 voltage wave V1-. At the customer premises CP, a communication terminal 
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(e.g., a phone set) is wired to the termination TERM 2 of the local loop LP. 
The communication terminal is modeled as a variable impedance ZL. 

A first reflectometry measurement takes place while the phone 
line is idle (see fig. 3). The phone set at the customer premises CP is on-hook 
5 and acts as an open circuit (ZL = oo). The prediction of the loop topology and 
the prediction of the network parameters of the local loop are carried out 
therefrom. 

The prediction of the loop topology yields the following results 
(refer to the above mentioned patent for further details): 
10 - a first line section {z oi9 e~ rih ) with a splice SP at 100 \xs from the 
termination TERM 1, 

- from the splice SP on, 2 line sections (z o2 ,e~ r2/2 ) and {z o3t ,e~ nh ) with 
respectively 75 \xs and 25 \xs as transmission delay. 

Once the tree structure of the local loop has been determined, the 
15 corresponding parametric model is selected from a predefined collection of 
parametric models. If no match is found, the operator is requested via the 
output device OUT to import or to create the corresponding parametric 
model. 

The 1-port scattering parameter S™ odel for a subscriber loop with a 
20 bridged tap is given by: 

smodel = -p.+p^™ + Pg p 2 p„e- 2 ^ +p t/W T»* -/W- 2(r,Wi) -Ka^^-p.pMp^^+pMe™™™ 
1 - P 9 Mf** ~ P2P„e- 2 ^ - P 3 p l2 e' 2 ^ + Pg p 2 p l2 e-™ + ^ + p gP ,p u e™^ + p 2 p,p 2 ,e^^ - p g p 2 p,e^r^ 

where: 

. P = z *- z " 

Ps z g + z 01 ' 

- P tj = l 0i l 0J ~l 0i l 0k ~l° jZ 7 0k , with (i,j,k) = {(1,2,3),(1,3,2),(2,3,1)}, 

25 - p, = §— with i = {2,3}, 
z i + z 02 

- Zi is the impedance load at the termination TERM i. 
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The transmission line type, as known to the telecommunication 
operator and as entered via the input device IN, is an Unshielded Twisted 
Pair (UTP) consisting of 2 parallel circular cylindrical conductors immersed 
in a homogeneous dielectric, which is assumed to be linear and isotropic. 
(z 0 ,e-*) of a UTP is given by: 
yi 



where: 

- s = j co = jlrf represents the Laplace frequency, 

10 - *F = 3a 2 J^J 2 + 2a 2 J x J 2 + a 2 J 0 J 3 ci eno tes an auxiliary function used in 

a 2 J 2 J 3 + tf 2 «/,./ 2 + 3a 2 J 0 J 3 + 

the calculations, 

- J, = Ji{a^4-s) denotes the Bessel function, 

- y p =a 5 s/l is the parallel per-unit length admittance of the transmission 
line. 

15 The coefficients ai to as may be given initial values in accordance 



with: 



1 jfi € t 



an\<7 arcosh(£)/2a) 



D 

a 3 = a^jjua 



20 a A = 



4 ln(l/V^") 



arcosh(D/2a) 



where: 
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- a is the radius of the circular conductors (same radius for both 
conductors), 

- D is the distance between the 2 conductor axis, 

- a is the conductivity of the conductor material (e.g., copper), 

5 - 8 and ^ are respectively the electrical permittivity and magnetical 
permeability of the dielectric material (e.g., polyethylene). 

15 coefficients are thus defined for the parametric model of the 
local loop LP, 5 for each line section. The coefficients need to be optimized 
with respect to the first measured reflectogram. Refer to the above 
10 mentioned patent for further details about the cost function and the 
optimization algorithm. 

Once the coefficients have been optimized, the theoretical channel 
capacity of the local loop can be computed from the formula: 



C=f log 3 



df 



15 where: 

- H(f) denotes the transfer function between the central office CO and the 
customer premises CP, 

- S(f) denotes the power spectral density of the signal transmitted at the 
central office CO over the local LP, 

20 - N(f) denotes the power spectral density of the noise at the customer 
premises CP. 

As expected, the following issue arises: which one of the transfer 
functions H12 or H13 is the transfer function between the central office CO 
and the customer premises CP or, alternatively, which one of the loop 
25 terminations TERM 2 or TERM 3 the customer premises CP are associated 
with? Therefore, a second reflectogram is executed to determine the loop 
termination associated with the customer premises CP. 

The second reflectometry measurement takes place while a phone 
call is on-going (see fig. 4), the phone set at the customer premises CP being 
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off-hook (ZL * co). The phone call can be initiated by the technician himself 
or by a third party. 

The first time domain reflectogram is then compared with the 
second time domain reflectogram in the processing unit PROC by 
5 subtracting one from the other, that is to say: 
A(n) = \s x (ri)-s 2 (n)\ 9 

where: 

- si(n) is the n th averaged sample of the first time domain reflectogram, 

- S2(n) is the n th averaged sample of the second time domain reflectogram. 
10 The result of the subtraction is smoothed by correlating A with a 

Gaussian function h: 

+00 

A\n)= £A(n-*)x /*(*:) 

Ar=-oo 

where: 

- *(*) = 3-7,(0, 

15 - I k (t) denotes the modified Bessel function of integer order k, 

- tis the scaling factor, to be chosen adequately. 

A' function is finally used to determine the round trip delay 
between the measurement side and the customer side. First, the global 
maximum M of A' function (dA'/dt = 0) is looked for (see fig. 5). Next, the 

20 first inflexion point I of A' function (d 2 A'/dt 2 = 0) preceding the global 
maximum M is looked for (see fig. 5). The time corresponding to this 
inflexion point with respect to the time at which the initial signal has been 
coupled to the local loop LP defines the reference delay DREF (see fig. 5). 

The reference delay DREF is used ultimately to identify the loop 

25 termination associated with the customer premises. This is achieved by 
selecting, out of all the paths between the loop termination at which the 
measurements are performed and any other loop termination, the path 
which the accumulated delay matches the closest the reference delay DREF. 
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Still referring to fig. 2, the round trip delay between TERM 1 and 
TERM 2 is (100 + 75)x2 = 350/tf and the round trip delay between TERM 1 
and TERM 3 is (100 + 25)x2 = 250/^ . By comparing them against DREF = 

350 (is, TERM 2 is undoubtedly the loop termination associated with the 
5 customer premises CP. 

Once the loop termination associated with the customer premises 
CP has been identified, the proper transfer function is selected for the 
computation of the theoretical channel capacity between the central office 
CO and the customer premises CP. 
10 The transfer function between the central office CO and the 

customer premises CP is thus given by: 



// = //„=- 



r 23 r / (g"^^ +r *' 2 } + g" (r|/ ' +ri ' 2 +2r *' 3 ) ) 



where: 

- r 9 = 1 + Aj . with (i j) = {(1,2),(1,3),(2,3)}, 



15 - A = 



Z -7 



The theoretical channel capacity between the central office CO 
and the customer premises CP is computed therefrom, with Zi and S(f) 
being substituted by their respective service value (e.g., Zi = 100 or 135 Q 

20 for a DSL modem). The so-computed capacity is eventually displayed via the 
output device OUT as a yes/no answer to whether the loop capacity is 
greater or less than a pre-determined target figure, that is to say as a yes/no 
answer to whether the local loop is qualified for that figure or not. 

In an alternative embodiment of the present invention, the input 

25 device IN could comprise any of a variety of input devices as known to the 
skilled person, such as dedicated push buttons, a touch screen, etc. 
Similarly, the output device OUT could comprise any of a variety of output 
devices as known to the skilled person, such as a printer, an output file, etc. 



- 16 - 



In an alternative embodiment of the present invention, the time 
domain reflectometry unit TDR would monitor the line status by measuring 
the power spectral density of the local loop in the voice band and by 
determining whether it is greater or less than a pre-determined threshold. 
5 The line status would be latched in a register that is addressable by the 
processing unit PROC via the I/O address range. The processing unit PROC 
would be interrupted with a dedicated interrupt vector whenever the line 
status changes. The processing unit PROC could also polls cyclically the line 
status register. The time domain reflectometry unit TDR could also notify 
10 the line status changes to the processing unit PROC via the data 
communication bus COM by means of dedicated message primitives. The 
processing unit PROC would automatically trigger a first time domain 
reflectometry measurement while the line is idle (PSD < threshold) and a 
second time domain reflectometry measurement while the line is busy (PSD 

15 > threshold). In so doing, the processing unit PROC would ascertain the line 
status during the reflectometry measurements and the whole qualification 
procedure would be automated. 

In an alternative embodiment of the present invention, the input, 
output and processing means would reside on a remote host, while the time 

20 domain reflectometry unit TDR would be mounted on a card and plugged 
into a card slot of a telecommunication network element such as Digital 
Subscriber Line Access Multiplexer (DSLAM). The time domain 
reflectometry unit TDR would be coupled to the processing means via a data 
communication network such as the internet, a Local Area Network (LAN) 

25 or a Wide Area Network (WAN). 

A final remark is that embodiments of the present invention are 
described above in terms of functional blocks. From the functional 
description of these blocks, given above, it will be apparent for a person 
skilled in the art of designing electronic devices how embodiments of these 

30 blocks can be manufactured with well-known electronic components. A 
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detailed architecture of the contents of the functional blocks hence is not 
given. 

While the principles of the invention have been described above in 
connection with specific apparatus, it is to be clearly understood that this 
5 description is made only by way of example and not as a limitation on the 
scope of the invention, as defined in the appended claims. 



